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Abstract. Within the framework of hierarchical cluster- 
, ing scenarios, we investigate the consequences for the 
properties of virialized halos of the constraints provided 
, by numerical simulations on the first few correlation func- 
tions. Thus, we show that the density field cannot be 
, described by a collection of smooth halos with a univer- 
■ sal density profile. This implies that substructures within 
larger objects play an important role (but a mean spher- 
ically averaged density profile may exist). In particular, 
a possible interpretation is that collapsed objects can be 
divided into an infinite hierarchy of smaller objects with 
increasingly large densities (these substructures might also 
be continuously destroyed and created by the long-range 
action of gravity). Finally, we present multifractal mod- 
els (restricted to non-linear scales) which can describe in 
a natural way such non-linear density fields with increas- 
ingly large fluctuations at smaller scales. We relate their 
\ properties to the correlation functions and present a few 
constraints they are expected to satisfy, using theoretical 
considerations as well as constraints from numerical sim- 
ulations. Thus, the simplest realistic model is the bifractal 
model described in Balian & Schaeffer (1989a). Moreover, 
we show that it should provide (at least) a very good ap- 
proximation of the multifractal properties of the actual 
non-linear density field, hence of the probability distribu- 
tion of the density contrast. The implications of this model 
(e.g. for galaxies) are detailed in other studies. 
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1. Introduction 

In the standard cosmological scenario the non-linear gravi- 
tational structures we observe in the present universe have 
formed by the amplification of small primordial density 
fluctuations. Moreover, according to hierarchical models 
the power increases at small scales so that low masses 



collapse first to build small virialized objects which will 
later merge to form broader halos as larger scales turn 
non-linear. The objects built by this continuous process 
produce the galaxies and clusters we observe today. Thus, 
the description of the non-linear evolution of the density 
field is an important task in cosmology. Most works have 
considered three tools to measure its properties. First, one 
can study the mass function of halos defined for instance 
by a given density threshold. This is the goal of the pop- 
ular Press-Schechter approximation (Press & Schechter 
1974). Note however that this analytical approach is re- 
stricted to only one value of the density contrast threshold 
(A ~ 177). Second, one can consider the density profile of 
individual objects. This is mainly done through numerical 
simulations (e.g. Navarro et al.1996; Tormen et al.1997). 
Third, one can measure the many-body correlation func- 
tions. These are closely related to the counts-in-cells: the 
universe (or the simulation box) is divided into cells of a 
given size and one considers the probability distribution 
of the density contrast realized in such a cell. This can be 
studied through numerical simulations (e.g. Colombi et 
al.1997; Valageas et al. 1999a) and analytical works based 
on specific scaling models (Balian & Schaeffer 1989a). 

While the first two properties are directly linked to the 
characteristics of the astrophysical objects we observe the 
third has the advantage to be better defined (it bypasses 
the problem of the recognition of individual halos) and to 
contain much more information. In particular, it should be 
possible to obtain the properties of individual halos from 
the correlation functions. For instance, the mass functions 
of objects defined by various density thresholds (not nec- 
essarily constant) obtained from the scaling model devel- 
opped by Balian & Schaeffer (1989a) for the correlation 
functions were considered in Valageas & Schaeffer (1997) 
and Valageas et al. (1999a). 

In this article, we investigate the general link between 
the behaviour of the correlation functions and the shape 
of collapsed halos. First, in Sect. 2 we introduce our nota- 
tions. We relate the correlation functions to the moments 
of the density probability distribution and we note a few 
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exact properties. In Sect. 3 we recall the constraints on 
the first few correlation functions provided by numerical 
simulations. Then, in Sect. 4 and Sect. 5 we show that the 
density field cannot be described as a collection of smooth 
halos with a universal profile. This points out the role of 
substructures, as discussed in Sect. 6. Finally, in Sect. 7 we 
present multifractal models which provide a natural tool 
to describe a density field with increasingly large density 
fluctuations at smaller scales. We discuss their link with 
correlation functions and we obtain some constraints from 
theoretical considerations as well as from numerical sim- 
ulations. In particular, we show that a bifractal model 
should provide (at least) a very good approximation for 
the non-linear density field. 

2. Scaling exponents 

A convenient description of the density field is provided by 
the probability distribution Pi(pi) dpi of the local mean 
overdensity over spheres of radius I and volume V\. Here 
we defined: 



p(r) = ^ = (l + S)(r) 
and 

Pl(r)= / -jrP'r' >0 
Jvi n 



(1) 



(2) 



Since usual probability distributions can be denned by 
their moments one is led to consider the parameters p P (l) 
given by: 



/>oo 

peH : p p (l) = (pf) = / dppPPdp) 

Jo 



(3) 



where () denotes an average over the position r of the 
sphere. Of course, if all moments p p (p <G HI) are finite 
their values for p E IN (or — p <E IN) are sufficient to define 
the density field. Hence all parameters p p (p € IR) can be 
obtained from the sole p p realized for p e IN (or — p € IN) . 
From the definition (3) one can easily show that: 



p < p > , q > : to < ^9 
Pp Ppi 

In particular, we obtain: 
p e JR : to > 



Pp Pp—l 

Since p\ = 1 we have for positive integers: 

p e IN* : f J "P+ 1 > „ 2 anc [ /ip > uS 1 
P P 



(4) 



(5) 



(6) 



We also obtain, using po = 1 (and the fact that (1 + 5) 
and (1 + <5) _1 are not proportional), 



P-(p-i) 



(7) 



hence: 

p G IN* : p- p > p v _ x 

From the definition (3) we can also write: 

peW : p p (l) = d " ri ^ p{v l )...p{vp) 



(8) 



(9) 



which emphasizes that these moments are related to the 
p— point correlation functions £ p . We define: 

lp{l)=l d \f rp ep(ri,-,r p ) , ? =?i = l (10) 

JV l Vl 

Then, from (9) and the definition of correlation functions 
(Peebles 1980) we have: 



PGIN* : p p (l)= J2 Oii.-.ip ^(0-^(0 



(11) 



where the sum runs over the integers < i\ < ... < i p 
with i\ + ... + i p = p, the coefficients a^^.^ p are positive 
and ai,..,i = ao,..,o,p = 1- In particular, 



Pi = l + £ 



(12) 



where we note £ 2 as £• As a consequence, in the non-linear 
regime we obtain from (4) and (11): 



p E IN , £ — > oo : 



Mp 



(13) 



Of course this also implies that the density field is strongly 
non-gaussian in the non-linear regime. We are interested 
in the case of a critical universe with an initial power- 
spectrum which is a power-law: P(k) oc k n . Then, it is 
natural to expect a power-law behaviour for the various 
quantities £ p (Z) and p p (l) in the non-linear regime £ 3> 1 
(i.e. at small scales). Note indeed that numerical simula- 
tions show a power-law behaviour for £ (Jain et al.1995; 
Colombi et al.1996; Valageas et al. 1999a). Thus, we define 
the limits: 



p G IR : 7 P = - hm p p v 
F i^o In/ 



(14) 



From (13) this means that for integers p we have: 

P G IN : £ p (Z) - r 7 " for / -» (15) 

where we ignore logarithmic corrections. From (5) and (7) 
we obtain: 



p G IR : 



7p+i - 7 P > 7 P - 7p-i 



(16) 



,7-(p+i) - 7-p > 7-p - 7-(p-i) 
In particular, we have: 

'7p > (P~ 1) 72 



7o = 7i = , p G IN* : 



(17) 



7-p > P7-1 
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Moreover, since p = (1 + 5) is a positive quantity l 3 pi(r) is 
a non-decreasing function of I for any point r. This implies 
that for any p > the average l 3v (pf) is a non-decreasing 
function of I, so we get: 

p > : 3p - 7 P > hence lim ^ < 3 (18) 

p^oo p 

Note that all the results obtained above are only due to 
the fact that p > 0. To obtain the probability distribution 
Pi(p) dp of the overdensity at scale I it is convenient to 
introduce the quantities: 



peJN* : S p (l) 



and to define the function: 



Si = 1 



vi{y) 



(19) 



(20) 



P =i 



Then, one can show (White 1979; Balian & Schaeffer 
1989a) that: 



Pi(p) 



/-\-ioo 
-zoo 



i^L e (py-vi(y))/6 



(21) 



The interest of the function (y) is that it shows a very 
small (if any) dependence on scale I in the non-linear 
regime, so that the evolution with I of the probability dis- 
tribution Pi(p) is given by the sole two-point correlation 
function £. This was shown in the non-linear regime from 
numerical simulations by Colombi et al.(1992, 1994, 1995, 
1997), Munshi ct al.(1999) or Valageas et al. (1999a) for 
instance. Note that this property also holds in the quasi- 
linear regime (Bcrnardeau 1994), with a different function 
(pi (y) . In the regime £ 3> 1 which we consider in this arti- 
cle, we have from (19) and (17) for p > 3: 



I -» 



s p (i) ~ r 



Cp = 7p-(p-1)72 >0 (22) 



Thus, the coefficients S p are constant or increase with £ 
(neglecting logarithmic corrections). We recall in appendix 
the behaviour of Pi(p) and the moments p p implied by 
simple forms of <pi(y). Note that one may also define from 
the moments \x p the generating function: 



oo (_-i\p 

My) = Y.—r^ l )y 



p=0 

which leads to 



and: 

My) 



-Mi y)/S 



(23) 



(24) 



(25) 



Developping (25) in y one can obtain the coefficients 
«»!,..,*„ defined in (11). 



3. Constraints 

The behaviour of the two-point correlation function £ and 
of the first coefficients S p with p < 6 has been studied 
by many numerical works in the non-linear regime. Thus, 
the slope 72 (usually noted 7) of £ is consistent with the 
prediction of the stable-clustering ansatz (Davis & Peebles 
1977; Peebles 1980): 



CI 



3(3 + n) 
5 + n 



where n is the slope of the initial power-spectrum (Jain et 
al.1995; Colombi et al.1996; Valageas et al.l999a). More- 
over, the stable-clustering assumption predicts that in the 
highly non-linear regime the p— point correlation functions 
satisfy the scaling law: 

e p (Ar 1; Ar p ; a) = a 3 ^ A"^" 1 ) | p (n, ...,r p ) (26) 

where a(t) is the scale-factor, which means that ( p — and 
the coefficients S p do not depend on the scale /. Colombi 



et al.(1996) found a very small drift S p oc £ 



- 0.045 (p-2) 



for 



p < 5 while Munshi et al.(1999) obtained no scale depen- 
dence for p < 6. Thus, numerical simulations provide the 
conservative constraint: 

C2 : < C P < 0-045 (p - 2) 7 for 3 < p < 5 

Hence, the first scaling exponents 7 P are very close to the 
stable clustering prediction 7 p c : 



3 <p< 5 : < 



lp-T P 



< 0.045 



p-1 



with p > 1 : 7p C = (p - 1) 7 



(27) 



(28) 



where 7 is given by CI. The results CI and C2 pro- 
vided by numerical simulations are certainly robust. In- 
deed, they arc derived (among other methods) from the 
counts in cells which are a very well defined procedure. 
By contrast, we think that the properties of dense halos 
seen in numerical simulations may be more uncertain since 
it is not always easy to identify these objects which may 
not be spherical. However, many numerical studies have 
considered the density profile of virialized halos. In partic- 
ular, Navarro et al.(1996, 1997) found that the spherically 
averaged density profile follows a universal behaviour, in- 
dependent of time and of the initial power-spectrum (al- 
though the numerical parameters which enter this average 
profile depend on the index n of the power-spectrum its 
shape does not depend on n). Thus, many studies (e.g. 
about galaxy formation processes) assume that: 

HI : the density field can be described as a collection 
of spherical halos with a smooth density profile given by a 
universal law. 
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We note that this simple (and convenient) picture ne- 
glects substructures within large halos. Thus, HI is more 
restrictive than the assumption of a spherically averaged 
universal profile which removes substructures by defini- 
tion. However, the validity (or non-validity) of HI will 
provide some information on the importance of substruc- 
tures. Wc note that numerical simulations often suffer 
from the overmerging problem: due to poor resolution 
small objects often get disrupted within large halos. Nev- 
ertheless, some simulations with a high resolution show 
significant substructures (Moore et al.1998; Ghigna et 
al.1998). In the next sections, we shall check whether the 
description HI is consistent with the constraints CI and 
C2. 

4. Single power-law profile 

We first consider the simplest case where the density field 
can be described as a collection of spherical halos with a 
smooth power-law density profile and a fixed overdensity 
(1 + A c ) at their radius R (e.g. A c <~ 200 from the spher- 
ical collapse model). We note the multiplicity function of 
virialized halos of mass M, volume V(M), as r](M)^f 
while the local overdensity at distance r from the center 
of a halo of radius R is: 



p h (r) = (1 + A C ) (L) 



-0 



(29) 



where < (3 < 3 is the slope of the density profile. The 
mass function satisfies: 



f°° M ,*^ dM , 



(30) 



which states that the fraction of matter enclosed in these 
halos is unity. It could be smaller if one adds for instance a 
uniform component. However, this would not change our 
results since they only rely on the fact that the integral 
in (30) is finite and this low density component only con- 
tributes by a constant amount to the moments fi p (p > 2) 
which tend to infinity when I — > 0. We want to obtain the 
exponents j p for p > 2 realized in such a density field, to 
compare them with the constraints CI and C2. Thus, for 
the sake of clarity we shall not consider the normalization 
factors which enter our calculations since we only need the 
exponents of the various power-laws of interest. To obtain 
the behaviour of p, p we must separate large halos (i? > I 
one cell is within one such halo) from small halos (i? < I 
one cell may contain several whole halos). Thus, we write 



/ip(0~Ai(0 + A 2 (0 



(31) 



where A\ is the contribution of cells which are enclosed 
within a massive halo (R > Z) while A 2 is the contribution 
of cells which contain zero, one or several small objects. 
Note that for R <~ I our distinction is not adequate but 
we will not consider in detail such intermediate regimes 



which do not influence the power-law behaviours we are 
interested in. Then, 

f 00 dM 

Ml)~J V(M) — V(M)(pP)i,M (32) 

where (p p )z,m is the mean overdensity to the power p seen 
by a cell I which intersects a halo M. Thus, 



M 



1 

V 



J\nr 2 drp h (rf + ^p h (lf 



which leads to: 
(p p )lm~ (1 + A c )p 



1+[ R 



3-p/T 



(33) 



(34) 



As explained above, here the factor 1 stands for a numer- 
ical constant. Thus, we have: 



/ V(M) 

JR=l 



dM M 
M (l + A c )/5 



(i + A c y 



3-P/3' 



Considering separately the cases (3 — p(3) < and (3 — 
p(3) > 0, and using (30), we can push the lower bound of 
the integral down to and we obtain: 



I -» 



A^l) ~ P + f-Pt 3 



(35) 



Here, we also assumed that J °° M 1+a Tj(M)dM/M con- 
verges for any a > 0. This means that the multiplic- 
ity function has an "exponential" cutoff at larges masses 
which is consistent with numerical results and observa- 
tions of the galaxy luminosity function (of course the ar- 
gument of the exponential is not necessarily linear in M). 
This ensures that p p (l) is finite for any p > 2 for a finite 
I > 0. On the other hand, the overdensity within cells 
which contribute to ^(Z) is smaller than (1 + A c ) (which 
is realized when the whole cell is covered by small halos) 
which means that A 2 (l) has a finite limit for I — > 0. Thus, 
we obtain: 

p>2 : p, p (l) - Z° + l 3 - pf) (36) 
If this density field is consistent with CI then: 

/»-»±2.»filM (37) 
2 n + 5 

as was obtained by McClelland & Silk (1977), and more- 
over 



P > 2 : 7 P = p(3 - 3 



C P = (P- 2) 



3-7 



(38) 



Thus, the exponents ( p violate the constraint C2 for any 
power-spectrum of interest — 3 < n < 1. This shows that 
the model considered in this section cannot provide an ac- 
curate description of the actual density field. Moreover, we 
note that the slope (3 needed to reproduce the two-point 
correlation function £ is much larger than the density pro- 
file observed in numerical simulations (/3 = 1 in Navarro 
et al.1997; (3 = 1.4 in Moore et al.1998). 
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5. Mass-dependent profile 

While the density profile of the halos considered in the 
previous section was given by a unique power-law for all 
objects, many authors have found that small objects are 
denser than larger ones. Thus, Navarro et al.(1996) obtain: 



Ph(r) = 



{r/R s ){l+r/R s f 



with 



S c = 



200 c 3 

~Y ln(l + c)-c/(l + c) 



(39) 



(40) 



where R s = R/c is a characteristic radius (c > 1) and 
small objects are more centrally concentrated: they have 
a higher c and 5 C . The same trend is seen (Tormen et 
al.1997) using a Hernquist profile (Hernquist 1990): 



Ph{r) 



(r/R a )(l + r/R a f 



(41) 



Thus, these objects have a mean overdensity (1 + 6) <~ 5 C 
over the radius R s and (1 + 5) ~ 200 over R. Navarro et 
al.(1996) have found that the characteristic overdensity of 
these halos is proportional to the mean density of the uni- 
verse at the epoch they were formed (which we note by the 
rcdshift zj): S c oc (1 + zj) 3 . This redshift Zf increases for 
smaller halos which formed when the corresponding scale 
M became non-linear. Indeed, using the extended Press- 
Schcchtcr model (Press & Schechter 1974) developped by 
Lacey & Cole (1993, 1994), one obtains good results by 
setting: 



erfc 



1.69 z f (M) 



y / 2(a 2 (M/2) — a 2 (M)) 



(42) 



where as usual a 2 (M) is the variance of the density fluc- 
tuations at scale M in the present universe given by the 
linear theory (Lacey & Cole 1994). Note however that 
Navarro et al. (1996, 1997) get better results by using a 
small mass fraction / = 0.01 rather than / = 1/2 in (42). 
For small masses, this leads to: 



M <C M* : (1 + z f ) 



„<M> - , £ 



-(n+3)/6 



(43) 



where we defined by cr(M*) = 1. Hence, we obtain: 

(44) 



/ M \ -(n+3)/2 

M <C M* : 6c ~ ( jjr 



As we noticed above, this simply states that the charac- 
teristic density of halos of mass M <C M* is given by 
the mean density of the universe at the time when this 
mass scale became non-linear. This behaviour is probably 
more robust than the detailed prescription (42). We note 



that this means there has only been a negligible evolution 
of these small objects: their density has not significantly 
changed since their formation. Note that Salvador-Sole et 
al.(1998) also recovered (44) using a more detailed model 
than (42) which takes into account violent mergers and 
slow accretion. Neglecting the logarithmic correction in 
(40) we also have: 



M{< R s )~p5 c R 3 



pc 3 



pR 3 ~ M(< J?) (45) 



This means that most (or a finite fraction) of the mass 
is enclosed within R s . This would be exact for any outer 



density profile (r > R s ) steeper than ph(r) oc 



like 



the Hernquist profile (41). This is also consistent with the 
results of Navarro et al.(1996) since simulations cannot 
distinguish between ph(r) oc r~ 3 and Ph(f) oc r -31 (for 
instance). Thus, we shall consider in this section a model 
where large halos have a single power-law density profile: 

M > M* : p h (r) = (1 + A c ) (£) (46) 
while small objects follow a double power-law profile: 



M < M* 



Rs 
~R 



M 



(n+3)/6 



and 



(47) 



/ x -(n+3)/2 / n. 

r<R s : p h (r) = (1 + A C ) (^) 



M 



R s <r<R : p h (r) = (1 + A c ) y w 
with: 



-(n+3)/2 



(*)■ 



(3 < 3 



> 3 



(48) 



Of course, this double power-law profile is a simplification 
of the smooth profiles (39) and (41). However, it captures 
their main features which are that a finite fraction of the 
mass is enclosed within a radius R s with an inner density 
profile with a slope (3 = 1. The outer slope $ will not enter 
our final results (it only ensures that R s contains a finite 
fraction of the mass and that M oc R 3 ). Since small halos 
may provide a significant contribution to the moments fi p 
(contrary to the case considered in the previous section) 
because of their higher density, we must detail the be- 
haviour of the multiplicity function in the low mass limit. 
Thus we write: 



M <C M* : r](M) oc M 



o-i 



9>0 



(49) 



where we did not consider logarithmic corrections. The 
lower bound 9 > is given by the condition (30) . In most 
models, one also has 9 < 1 so that the number of small 
objects diverges (as in the Press-Schechter mass function, 
Press & Schechter 1974). We note that this is also nec- 
essary to reproduce the faint-end slope of the galaxy lu- 
minosity function, since the luminosity usually scales as 
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a power of the dark matter halo mass in order to ob- 
tain the Tully-Fisher relation (although the actual power- 
spectrum is probably not a power-law it is smooth enough 
to be approximated by a power-law over the mass range 
corresponding to galaxies). However, we shall not need 
this upper bound in the following. The fraction of mat- 
ter (or luminosity) contributed by small objects in these 
models and observations also decreases for small M (or 
L) as a power- law with a positive exponent, that is 9 > 0. 
Thus wc do not consider here the case 8 — which can 
still be made consistent with (30) with suitable logarith- 
mic prefactors (e.g. rj(M) oc l/[Mln 2 (M)]). Hence, this 
model corresponds to the hypothesis HI where the profile 
of the halos is given by the results of Navarro et al.(1996) 
or similar studies (Hernquist 1990; Moore et al.1998). In 
a fashion similar to the calculations shown in the previous 
section we can obtain the scale dependence of the moments 
\i v for p > 2. Thus, we write: 



n„(i) - Mi) + Mi) + Mi) + Mi) 



(50) 



where A\ corresponds to cells embedded within a large 
halo M > M*, A2 to cells enclosed within an object M < 

such that I < R s , A 3 to cells within a halo M < 
with R s < I < R and A4 to cells which contain zero, one or 
several small objects. Indeed, since we are interested in the 
small scale behaviour we have I <C i?* (where i?* defined 
by a(i?*) = 1 is of the order of the radius of collapsed 
halos of mass M„). We obtain: 

+ , A 2 ^z° + ; 3 - p/3 + ^ [e - (p - 1)(ri+3)/21 

A 3 ~ zt&[9-(p-i)(»+3)/2] , At finite 
Hence, the exponents 7 P are given by: 

3(3 + n) 



P>2: 7 P = Max 



p/?-3,(p-l)- 



66» 



5 + n 5 + n 



(51) 



Since 9 > 0, the constraint CI implies (37) and (38). Thus, 
we recover the results of the previous section. This means 
that this model cannot provide a satisfactory description 
of the density field either. 

6. Substructures 

The results of Sect. 4 and Sect. 5 show that the assump- 
tion HI is not valid: one cannot describe the dark mat- 
ter density field by a collection of halos with a smooth 
density profile. In other words, the substructures embed- 
ded within larger halos play an important role: they have 
to be taken into account in order to obtain the correct 
scaling exponents 7 P . This is also directly seen from (51). 
Indeed, we note that both constraints CI and C2 are sat- 
isfied if 9 = 0. As we noticed in the previous section, this 
cannot be used as a satisfactory mass function of distinct 
halos (even after addition of ad-hoc logarithmic factors 



to ensure convergence) because this would contradict the 
observed galaxy luminosity function and the mass func- 
tions measured in numerical simulations. However, the 
value 9 = has a very simple and natural interpretation: 
one counts the same matter (particles) at small scales as 
is counted at larger scales. In other words there is some 
multiple counting and the same mass is seen alternatively 
as objects of scale M and smaller "sub-objects" of scale 
M 1 <C M. Indeed, let us consider a discrete model where 
a mass M (per unit volume) is recognized as A/i objects 
of mass Mi = X l M with A < 1 (e.g. A = 1/2) at any scale 
i > 1. In the logarithmic interval dM/M = A In M = In A 
there are A/j = M/Mi objects hence r?(M) oc 1/M and 
9 = 0. Thus, the model described in the previous section 
can be made to satisfy CI and C2 if we give the mass 
function i](M)dM/M a new interpretation: it now counts 
the substructures of mass M embedded within larger halos 
as well as new isolated objects. Moreover, the latter pro- 
vide a negligible contribution to the moments \i v (p > 2) 
as compared to the former. Hence, the previous consider- 
ations lead to a new simple and natural model: 

H2 : dark matter halos can be divided into substruc- 
tures of arbitrarily small mass M with a characteristic 
density equal to the mean density of the universe at the 
time when this mass scale became non-linear. Moreover, 
in this case the exponents j p (p > 2) are given by the 
stable- clustering approximation (28). 

Of course, logarithmic factors may be present: the 
mass embedded within substructures of scale lower than 
M may slowly decrease as l/|lnM| for instance. One 
must note that H2 docs not necessarily contradict the 
result from numerical simulations that the spherically av- 
eraged density profile of dark matter halos follows a uni- 
versal shape. Indeed, in such measures one automatically 
discards substructures in order to obtain a "mean" be- 
haviour. However, it is surprising (and somewhat suspi- 
cious in our view) that this "average profile" does not 
depend on the slope n of the initial power-spectrum while 
(at least) the first few correlation functions (p < 6) do de- 
pend on n (indeed, as emphasized by Navarro et al.1997 al- 
though the characteristic density S c of M* halos decreases 
with n they find that dark matter halos obey the same 
profile (39) for all power-spectra). In particular, one would 
rather expect a density profile of the form p(r) oc r^ 1 . We 
can also note that there is still a significant uncertainty 
on the shape of these halos (Moore et al.1998). The de- 
scription H2 predicts a significant amount of substructure 
(even with a logarithmic correction) while simulations of- 
ten only have a moderate amount of "sub-objects" . How- 
ever, this overmerging declines with higher resolution so 
that it may be mostly a numerical problem due to low 
resolution (Ghigna et al.1998). In addition, as we noticed 
earlier the measure of the moments fi p is a better defined 
operation than the identification of halos and substruc- 
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tures so that the conditions CI and C2 are probably more 
robust than the results obtained for the characteristics 
of individual objects. Hence we think it is safer to rely 
on these two constraints than on observations about the 
behaviour of halos themselves. Thus, H2 appears to be 
the simplest description which is consistent with numeri- 
cal results, while HI is ruled out. We note that a simple 
model which leads to H2 was described in Valageas (1998) 
(this article also described the quasi-linear regime for the 
density and velocity fields). A similar picture of the den- 
sity field was also obtained by Balian & Schacffcr (1989a) 
starting directly from the assumption that all exponents 
7 P (p G IN*) satisfy the stable-clustering prediction (28) 
(see also Valageas & Schaeffer 1997). We must stress that 
although the most simple interpretation of H2 is that 
small clumps embedded within larger halos do not evolve 
much and keep roughly their initial density it can also de- 
scribe a much more dynamical picture. Indeed, individual 
clumps may continuously form and disappear through the 
exchange of particles while the statistical properties of the 
density field remain unchanged. 

7. Multifractals 

7.1. General model 

We have seen in Sect. 4 and Sect. 5 that one cannot describe 
the non-linear density field by a collection of smooth halos. 
Hence, we must turn to another description. The impor- 
tance of the scaling exponents 7 P and the conditions CI 
and C2 obviously lead one to consider multifractal mod- 
els (e.g. Balian & Schaeffer 1989b; Frisch 1995) to describe 
the non-linear regime. We first introduce the cumulative 
probability Pi(p) = J-^ Pi(p)dp (for instance). Then, for 
any scaling exponent a G IR we define the scaling dimen- 
sion F(a) < 3 (which can be negative and take the value 
-co) by: 



lim^2 = 3-F (a) 
i^o Inl 

In particular, we assume that: 



(52) 



H3 : the density field is multifractal on small non- 
linear scales I <C Z„ where: 



: Pi(l a ) ~ l 3 -^ 



(53) 



where I* defined by cr(Z*) = 1 corresponds to the 
transition between the linear and non-linear regimes. As 
usual we did not write possible logarithmic factors. Note 
that this definition of multifractality is less restrictive 
than the (more traditional) assumption: pi(r) ~ l a for 
r e V a C IR 3 with dimX> Q = F(a), for I -> 0. The scal- 
ing exponents 7 p defined by (14) are obtained from H3. 
Indeed: 

MO = / P P Pi((>)dp ~ J P a+3 - F ^ P*(a)da (54) 



where P*(a) gives the weight of the various exponents. 
Then, using the steepest descent method, we obtain: 



7 P 



min [pa + 3 — F(a)] = max [F(a) — 3 — pa] (55) 



Thus, 7 P (seen as a function of p e IR) and F(a) are 
related by a Legendre transformation. We note a p the 
scaling exponent which corresponds to j p . Then, we have 
from (55): 



P > P ■ ay < ol v 
We also get: 
p e M : F'(a p ) =p 
Using (18) we obtain: 

a > -3 



if F'(a p ) exists 



(56) 



(57) 



(58) 



In other words F(a) = — oo for a < —3. Since 70 = there 
exists at least one point a such that F(ao) — 3. Using 
71 = we also obtain: F(a) < 3 + a and this equality is 
realized at least for one point a\. Thus we have: 



F(a ) = 3 , F(ai) = 3- 



a\ 



Oil < < ao 



(59) 



We show in Fig.l an example of F(a) and how one can 
obtain the exponents j p from this curve. 



iF(a) 



3+a, 



-3 / a. 



3+Tp 



p a 



Fig. 1. Geometrical construction of the scaling exponents 7 P 
from the fractal dimensions F(a). 



The results of Sect. 4, where the density field was 
modelled as a collection of halos with a single power- 
law density profile, can be obtained from this multifrac- 
tal description. Indeed, using the usual definition, in this 
case there are only two fractal dimensions: a = with 
F(0) = 3 (which corresponds to points in the halo of the 
objects or within a uniform component) and a = —f3 with 
F{—[3) — (centers of the halos). Here we considered a fi- 
nite number of halos since small objects do not contribute 
to \x p for p > 1. Then one recovers (38) from (55). Note 
that in this case ao = a-i — 0. We have already seen in 
Sect. 4 that such a model is ruled out. 
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7.2. Bifractal 

The simplest model which is consistent with CI and C2 
is the case of a bifractal, where there are only two scaling 
exponents: a — u\ < and a = a$ > 0. Then, CI implies 
ot\ = —7 and we note a = 7^/(1 — u) with < u < 1. 
Thus: 



F(- 7 ) = 3-7 , F 

F(a) = —00 otherwise 
and: 

'p>l-u> : 7 P = (p - 1) 7 
^p < 1 - uj : Ip^-P^, 



7W 

1 - LU 



(60) 



(61) 



This bifractal model is consistent with the stable- 
clustering assumption (26), see (28). Moreover, it satu- 
rates the inequalities (17). We note that, starting from 
the behaviour (26) for the many-body correlation func- 
tions, Balian & Schaeffer (1989b) obtained the fractal di- 
mensions implied by simple models for the density field 
which are consistent with these scaling laws and they re- 
covered this bifractal model. More precisely, they found 
F(a) = —00 for a < —7 and a > 70;/ (1 — u), while in the 
intermediate range —7 < a < 7w/(l — u>) the dimension 
F(a) must be lower than (or equal to) the line joining 
(—7, 3 — 7) to (70;/ (1 — co), 3). Thus, these points play no 
role for the exponents j p which are still given by (61). 

On the other hand, we can note that a larger class of 
models (than the bifractal) is consistent with the stable- 
clustering prediction (28) for the correlation functions 
since one only needs: 

F(a) = -00 for a < -7 and F(- 7 ) = 3-7 (62) 

The exponents 7 P for small p are related to the behaviour 
for large y of the function <pi(y) introduced in (20). Thus, 
if we have: 



y 



00 



ip(y) ~ ay w with < lu < 1 



(63) 



where ip(y) = lim/.+o fi (y) we recover the second line of 
the system (61), as shown in appendix in (A12). On the 
other hand, if tp(y) <~ In 2 y for y — > 00 one obtains j p = 00 
for p < 0. This corresponds to density fields where there is 
a large probability to find very underdense regions while 
(63) leads to an exponential cutoff for p — > 0, see (A7). 

We note that this model is consistent with a universal 
spherically averaged density profile for virializcd halos of 
the form described in Sect. 4 or Sect. 5. It simply adds to 
the two scaling exponents —7 and 7w/(l — uj) the points 
{-/3,0) (for ScctTwith now < 7) or (-7,3-7-6^/(5 + 
nj) (for Sect. 5, but in this case this point does not show 
on the curve F(a) since 6 > 0). However, even if such a 



"universal shape" exists, it does not govern the exponents 
7 P which are given by the internal density fluctuations 
whose contribution dominates the moments [i v and the 
relations (61) are still valid. We also note that the scaling 
dimension F{— 7) = 3 — 7 which characterizes the matter 
distribution can lead to the filaments and walls seen in 
numerical simulations (e.g. Bond et al.1996) since 1 < 
3-7 < 2. 

7.3. Specific constraints 

In order to obtain more detailed characteristics of the 
curve F(a) we must explicitely consider the specific case 
we are interested in: the growth through gravitational in- 
stability of small gaussian density fluctuations, in a crit- 
ical universe. Although one should obtain the exponents 
7 P and F(a) from the Liouville equation, which is not an 
easy task, it is possible to derive a few constraints from 
simple and reasonable arguments. 

We first consider the behaviour of underdense regions, 
which governs the exponents 7 P for low p and F(a) for 
large a. We assume that strong underdensities ("voids") 
are created by the amplification of small initial negative 
density fluctuations. In particular, at a given scale R very 
low density regions ((1 + 5) — ► 0) come from the expansion 
of rare initial underdensities. In the early linear universe, 
density fluctuations at the mass scale M are characterized 
by the parameter v. 



St 



a(M) 



(64) 



where 5l is the linear density contrast. As seen in Valageas 
(1998) rare underdensities v <C — 1 have a spherically sym- 
metric profile. Moreover they tend to become increasingly 
spherical as they expand (Bertschinger 1985), hence we 
can use the spherical model to obtain their dynamics (see 
also Bernardeau 1994b) before their expansion is slowed 
down when they join with neighbouring underdensities. 
Thus, we write that the probability to find an underden- 
sity (1 + 5) — > in a cell of radius R is smaller than the 
probability we would obtain if all initial underdensities 
would follow the spherical dynamics expansion. Thus, we 
write: 



p " l<1 + A)i L irrmT) p ^- )di 



(65) 



where Pl(Sl)o18l is the probability distribution of the ini- 
tial linear density fluctuations: 



Pl(5 l ) 



1 



-Sl/(2a\R m )) 



2na 



(66) 



Here, R m plays the role of a mass scale so that R m = 
(1 + A) 1 / 3 i?. In the limit (1 + S) -> 0, the spherical model 
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relates the actual density contrast to the linear theory 
prediction by (Peebles 1980; Valageas 1998): 



-co 



(1 + *) 



20 
27 



-3/2 



(67) 



Thus, we obtain: 



Pr(< l + S)< 



-27/20 (l+<5) 



-2/3 



dSr 



27TCT 



20 A 

17 Ji 



3/2 



(68) 



x exp 



Si 



2<J 2 (R m ) 



which leads to: 



for p 



n+3 



> 1 : 



25 «+5 / I 

P i(< P) <^TJ^P~~ 



27V2tt 



x exp 



(69) 



n+3 



This means that limj_>o hiP;(< ^ a )/lnZ = oo in the case 
-a(l - n)/3 + (n + 3) < 0. From the definition (52) we 
get: 



F(a) 



-oo for a > 



3(3 + n) 
1 -n 



(70) 



In the framework of the bifractal model (60) it implies: 
5 + n 



fractal 



u < 



(71) 



However, it is possible to obtain another upper limit for 
u>i without assuming the density field to be a multifractal. 
Indeed, let us assume as in (A3) that ifii(y) is a power- 
law for large y with an exponent u>i which may depend 
on the scale I (the bifractal model is the peculiar case 
where and a/ are scale- invariant). Then the low-density 
cutoff of the density probability distribution Pi (p) is given 
by (A7). The constraint (69) applied in the limit p — ► 
implies (1 — oji)/lui > (1 — n)/3 hence we obtain the scale- 
independent upper bound: 



I < I* 



LUl < 



(72) 



Here oji is not necessarily related to a multifractal dimen- 
sion, it only characterizes the density probability distribu- 
tion seen at scale I. The existence of this upper bound is 
due to the fact that the probability to have very low den- 
sities decreases strongly as p — > 0. Indeed, in the linear 
regime the density field is very close to uniform (because 
of the sharp cutoff of the gaussian) and the amplifica- 
tion of underdensities is only a power-law of time in both 



the linear and non-linear regimes (but with different ex- 
ponents). Hence at late times the density field still "re- 
members" the initial gaussian cutoff. Our result is only 
based on the assumption that the most extreme "voids" 
come from the amplification of the most negative primor- 
dial density fluctuations which seems quite reasonable. We 
did not take into account shell-crossing since this would 
rather affect collapsed regions (it does not occur for the 
spherically symmetric voids considered here, sec Valageas 
1998). Note that for the values of n of interest the bound 
(72) is more constraining than (71). We compare in Tab.l 
the values of oj obtained in numerical simulations with the 
upper bounds (71) and (72). 



Table 1. Exponent lo for various indexes n of the 
power-spectrum. We compare the upper bounds (71) and (72) 
with numerical results from Colombi et al.(1997) (line C), 
Munshi et al.(1999) (line M) and Valageas et al. (1999a) (line 
V). 



n 


-2 


-1 





1 


(5 + ti)/6 


0.5 


0.66 


0.83 


1 


3/(4 - n) 


0.5 


0.6 


0.75 


1 


C 


0.3 


0.5 


0.65 


0.65 


M 


0.33 


0.4 


0.55 


0.7 


V 


0.3 


0.4 


0.45 





We can see that the numerical results are consistent 
with the upper bounds (71) and (72). More generally, al- 
though our argument was not completely rigorous the re- 
sults (70) and (72) should be quite robust. Indeed, one 
expects deviations from spherical dynamics to slow down 
the decrease of the "overdensity" (1 + 6) which leads to 
smaller a and lo. Indeed, a lower a or lo means less extreme 
underdensities. This may also explain why the values of to 
obtained from numerical simulations are smaller than the 
upper bound (72). 

Now, we show that the existence (for n < 1) of an 
upper limit smaller than unity for loi implies that loi is 
scale-independent within the framework of a multifractal 
model. Since a admits a finite upper bound the sequence 
ct- v has a finite limit a + for p — ► +oo, see (56). Moreover, 
we have from (55): 



(73) 



T-(p+i) - 7-p = a -{ P +i) + P ( a -(p+i) - a -v) 

+ F(a-( P +i)) -F(®-p) 

If F(a + ) is finite and a + is isolated, or F'(a + ) is finite, or 
F'(a + ) = +oo, there exists a real p + such that for any p > 
p + we have a_ p = a+. If F(a) is twice differentiable on 
an interval [d,a+[ with a < a + and F'(a + ) = — oo then 
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ot-p < a + for any finite p and using a Taylor expansion 
with remainder and (57) we obtain (7-( p +i) — 1-p) = 
a_( p+ i) + R p with < (a_( p+ i) — ct- p ). Thus, in all 
these cases we get: 



lim (7. 



(p+1) 1-p 

) = a + which is finite 



Using (Ai2) we have: 

M-( P +i)(0 



p — > +00 : In 



(i- P (l) 



l-u>l 



\np 



(74) 



(75) 



From the definition (14) of the exponents j p we also have: 



I < I* , p e M : In 



M-( P +i)(0 



H- P {1) 



In 



f-l—p (I* ) 



- (7-CH-i) - 7-p) In (£) 



(76) 



Here ln(/x p (Z*)) is the value obtained in I — I* by extrapo- 
lation from the non-linear regime I — > 0. Thus, we obtain: 



Z < Z* : 



1 - 1 - wj. 



■ Inl^ lim -y-(H-D-7-P 
J*/ P^°o mp 



(77) 



Here again is defined by the asymptotic non- linear 
regime. Using (74) we see that ui is constant, apart from 
a possible sub-logarithmic I— dependence (e.g. ln[ln(l/Z)]) 
which can be safely neglected for any practical purposes. 
Thus, the assumption of a multifractal together with the 
existence (for n < 1) of an upper bound for the scaling ex- 
ponents a is sufficient to show that the exponent uu of the 
power-law regime of <pi(y) at large y is scale- independent. 
However, the normalization factor o; may vary with I. In- 
deed, from (A12) and (14) we obtain: 



a*i 
ai. 



(l-o!)(a + -a ) 



_-(l-a,)(a4 



-<*o)/7 



(78) 



Thus, the measure of the normalization factor a ; allows 
one to get an important constraint on the curve F(a): it 
provides the value of the maximum scaling exponent a+. 
Indeed, we have: 



a > a + : F(a) = —00 



(79) 



Note that a + > a and a + is lower than the upper bound 
(70) while the upper limit (72) is irrelevant. Thus, the co- 
efficient ai remains constant (neglecting logarithmic cor- 
rections) or tends to as £ — > 00 (which means that un- 
derdense regions have an even smaller density). We note 
that numerical simulations (Colombi et al.1997; Munshi 
et al.1999; Valageas et al. 1999a) for the counts- in-cells are 
consistent with a constant value of a/ while the behaviour 



of the mass functions studied in Valageas et al. (1999a) 
might suggest a very slow decrease of a/ at smaller scales. 
From the dispersion of the measure of a; in these numeri- 
cal results and the small decline of the mass functions we 
can obtain an upper bound for a + using (78). We show 
our results in Tab. 2. (to get a+ we used the mean u> of 
the various simulations in (78)). 



Table 2. Exponents a for various indexes n of the 
power-spectrum. We compare the upper bound (70) with the 
mean value a given by numerical results from Colombi et 
al.(1997), Munshi et al.(1999) and Valageas et al.(1999a). We 
also present the maximum exponent a+ allowed by simula- 
tions, using (78). The exponent 7 is given by CI in agreement 
with numerical results. The upper bound for a_ is obtained 
from the measured dispersion of y s using (85). 



n 


-2 


-1 





1 


3(3 + n)/(l-n) 


1 


3 


9 


00 


CtQ 


0.45 


1.13 


2.2 


4.15 


a+ 


0.59 


1.39 


2.6 


4.42 


1 


1 


1.5 


1.8 


2 




1.15 


1.70 


2.07 


2.19 



Thus, the behaviour of F(a) in the domain a > olq is 
very close to the bifractal model. Indeed, numerical sim- 
ulations show that F(a) = —00 for scaling exponents a 
very close to ao- Moreover, even if there exist scaling ex- 
ponents slightly larger than ao (which may not be the 
case) they do not affect the slope w of the density prob- 
ability distribution while the normalization factor ai only 
shows a very weak (if any) scale dependence. We note 
that the value of a + obtained from numerical simulations 
through (78) is significantly lower than the theoretical up- 
per bound (70) and very close to a . This suggests that 
a + may in fact be equal to a . Indeed, if one measures di- 
rectly the curve F(a) in a numerical simulation, using (52) 
for instance, one would not be able to obtain F(ao) = 3 
and F(a) = —00 for all larger scaling exponents. Due 
to finite resolution effects there would be some disper- 
sion which would create some spurious exponents slightly 
larger than a - Moreover, it is obvious from (52) that the 
sharp F{a) which appears for instance in the bifractal 
model is only realized in the limit I — ■> 0: for any simula- 
tion where the available range of / is necessarily finite the 
ratio — \nPi(lmi n )/ In l m in will show a smooth cutoff. 

We now turn to high density regions (i.e. small a). Let 
us assume that after collapse mass condensations see their 



overdensity evolve as: (1 + 5) oc 



,3/3 



If these halos keep 



the same density (or radius) with time (i.e. no evolution) 
while they get embedded within larger objects we have 
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(3=1. Since Sl oc a (and collapse occurs for Sl ~ 1) we 
obtain for high densities: (1 + 5) <~ 5 L ^. Hence we get: 



(1 + *) 



— - 



i+3 _ n+3 



and: 



F(a) 



-oo 



for 



a < 



3(3 + n) 
~n + 3 + 2/[3 



(80) 



(81) 



Using the constraint CI we get /? > 1. Thus, there are 
two possibilities within this description: i) after virial- 
ization collapsed halos do not see their properties (den- 
sity, radius) evolve much (i.e. at most through logarith- 
mic factors) and the exponents 7 P (p > 2) are given by 
the stable-clustering assumption (28) or ii) some objects 
see their density increase with time as a power-law and 
the exponents j p (p > 2) are larger than the values (28). 
As larger scales become non-linear and small objects get 
embedded within increasingly massive halos, interactions 
between clouds (e.g. collisions) can affect their proper- 
ties. Although one might have expected these processes 
to decrease the mean density of small objects (through 
disruptions, relaxation towards the lower density of the 
larger host halo) the constraints CI and C2 show on the 
contrary that their mean density is constant (except for a 
possible logarithmic decline) or increases. This description 
assumes that one can follow the behaviour of individual 
mass condensations. In fact, this may not be the case as 
density fluctuations could be continuously created and de- 
stroyed, so that particles belong alternatively to overden- 
sities and underdensities of evolving magnitude. Thus, the 
behaviour of mass condensations could be more intricate 
than those of "voids" . In particular, one should explicitely 
take into account shell-crossing. However, we can directly 
obtain constraints on the curve F(a) from numerical re- 
sults. The location y s j < of the singularity of <pi(y), see 
(A4) and (20), is given by: 



Vs,l 



lim E S p^ 

P^oo S p+1 (l) 



(82) 



As explained in appendix we assume here that y s j is finite 
and non-zero. Then, we obtain from (22): 



y s (l) = y s (l*) lim ( — 



■Tp-7 



(83) 



where again y s (l*) is defined from the non-linear regime. 
Since a admits a finite lower bound (58) (only due to the 
positivity and additivity of the mass) the sequence a p has 
a finite limit — a_ for p — > +oo, see (56). From (55) we 
obtain: 



lim (7 p+ i - 7 P ) = a_ , 7 < a_ < 3 

p— >+oo 

in a fashion similar to (74). Hence we have: 



Veil) = y s (i*) 7- 



(84) 



(85) 



Thus, if Pi(p) shows a pure exponential cutoff at some 
scale I (i.e. y s (l) is finite and non-zero) it will display a sim- 
ilar exponential cutoff throughout the non-linear regime. 
Moreover, the strength y 8 (l) of this falloff declines or re- 
mains constant at smaller scales: the importance of ex- 
treme positive density fluctuations grows or is station- 
ary as one probes deeper into the non-linear regime (as 
compared with the predictions of the stable-clustering 
ansatz). From the dispersion of the measure of y s j ob- 
tained in numerical simulations (Colombi et al.1997; Mun- 
shi et al.1999; Valageas et al. 1999a) we obtain using (85) 
an upper bound for a_. We show our results in Tab. 2. 
Note that we have: 



a < —a- 



F(a) = -00 



(86) 



Thus, as was the case for the large a domain discussed 
above numerical simulations constrain the curve F(a) to 
be rather close to the bifractal model for a < —7. The 
possible scaling exponents do not extend down to a = —3 
and their minimum value — a_ is close to —7. This sug- 
gests again that the existence of scaling exponents below 
—7 may be due to finite effects (note moreover that the 
values displayed in Tab. 2 are upper bounds). Although we 
assumed y s j to be finite and non-zero the constraint on 
a_ shown in Tab. 2 also applies if the cutoff is not a "pure" 
exponential. Thus, if we write more generally: 



p>£ : Pi{~p)- _ 2 
with ki > we obtain: 



~ \ LU~ I — 1 

a s ,i P 



exp 



P 



x s ,i £ 



00 : p p <~ £ — - X 



(87) 



(88) 



The case Hi > 1 corresponds to y s j = —00 (no singularity) 
and Ki < 1 to y s ^ = 0. In a fashion similar to (77), using 
the constraint a > —3, we obtain that Hi is constant (ex- 
cept for sub-logarithmic terms). From (88) we also obtain: 



£ p^oo \ p J Mp 



(89) 



which is similar to (82). Since K; exhibits at most a sub- 
logarithmic dependence this implies: 



x s (l) = x s (h) — 



-("--7) 



(90) 



which is exactly (85) we got in the peculiar case k = 1 
since x s ^ = —l/y s j. Note that (90) does not depend on n. 
Moreover, numerical results constrain the argument of the 
exponential to be close to a linear function of the density: 
k ~ 1. 

Thus, we have shown that the curve F(a) is very sim- 
ilar to the bifractal model for a < ct\ and a > a . The 
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intermediate regime is more difficult to constrain since it 
depends on the details of the dynamics. For instance, one 
could add a point (0, 3) if a finite fraction of the volume 
is filled by a smooth component. However, such modifica- 
tions only influence the exponents 7 P for < p < 1, which 
corresponds to intermediate density fluctuations between 
the characteristic density of mass condensations (given by 

x Sj ;£) and "voids" (p <~ a]^ 1 ^ £ u ). One needs to 
solve the equations of the dynamics in order to get more 
precise information. However, it is clear from the previous 
results that the bifractal model provides a reasonable de- 
scription of the density field. Moreover, it is the simplest 
realistic model. 

We note that in the bifractal model the density prob- 
ability distribution Pi{p) in the very low-density regime 
(A7) also defines its behaviour in the intermediate regime 
(A6), and conversely, through the parameters a and u>. 
This also holds for more general models where <fi(y) — 
o-l y 1- " for large y. Indeed, as we noticed in Sect. 2 the 
coefficients S p (or p p ) for p e IN are sufficient to define 
entirely Pi(p) although they are mainly sensitive to posi- 
tive density fluctuations. However, the parameters to and 
a depend in a somewhat intricate fashion on the moments 
Hp or S p (p <E IN) through the asymptotic behaviour of the 
function ipi (y) while they can very easily be obtained the- 
oretically from the moments /i p {—p € IN) through (A12). 
On the other hand, the parameters k, lu s or x s j are di- 
rectly given by the moments fi p or S p (p € IN) through 
(88). Thus, both sequences have their own advantages. 

8. Conclusion 

In this article we have investigated the consequences for 
the shape of virialized halos of the constraints provided 
by numerical simulations on the first few correlation func- 
tions. Thus, we have shown that the density field cannot be 
described by a collection of smooth halos with a universal 
density profile. This does not imply that a mean spher- 
ically averaged density profile does not exist. It means 
that substructures within larger objects play an impor- 
tant role. In particular, we have shown that a possible 
interpretation of the constraints provided by numerical 
simulations is that dark matter halos can be divided into 
an infinite hierarchy of smaller objects. Then, the char- 
acteristic density of substructures of a given mass is the 
density of the universe at the time when this mass scale 
turned non-linear. However, these small mass condensa- 
tions may not be permanent entities as they could be con- 
tinuously created and destroyed through the long-range 
action of gravity. Finally, we have presented multifractal 
models which can reproduce the observed behaviour of 
the first few correlation functions. Thus, it appears that 
the simplest realistic model of the non-linear density field 
is the bifractal model developped by Balian & Schaeffer 
(1989a, 1989b). Moreover, we have shown that numerical 
results constrain the multifractal properties of the actual 



non-linear density field, and the probability distribution 
of the density contrast, to be very close (or identical) to 
the characteristics of such a bifractal. In the main text we 
have only considered the case of a critical universe with 
an initial power-spectrum which is a power-law. Never- 
theless, our main results (in particular the importance of 
substructures and small-scale density fluctuations) should 
also apply to power-spectra which are not pure power-laws 
but smooth enough to be reasonably approximated by a 
power-law over some significant range of mass (e.g. CDM) . 
However, in such a case the scaling exponents j p may ex- 
hibit a slow scale-dependence as a function of the local 
slope of the power-spectrum. Our results should also ap- 
ply to a low-density universe. In fact, since in such a model 
non-linear scales collapsed when the universe was still 
close to critical the non-linear regime should be the same 
as for the case ft = 1 with the same power-spectrum. One 
only needs to take care of the different time-dependence of 
the normalization of £ in the non-linear part when fl gets 
small, as done in Peacock & Dodds (1996) and Valageas 
& Schaeffer (1997). The implications for astrophysical ob- 
jects of this description have already been presented in 
other studies (e.g. Valageas & Schaeffer 1998 for galaxies; 
Valageas et al. 1999b for Lyman-a clouds). 



APPENDIX 

A. Density probability distribution 

Here we briefly recall the behaviour of the probability dis- 
tribution of the overdensity pi in spheres of size I implied 
by simple forms of (fii(y) (see Balian & Schaeffer 1989a for 
details). The relation (21) can be inverted as: 

_ z* 00 _ 
e -<Pi(y)/i = / e -~py/t p t (p)dp (Al) 
Jo 

This implies (using \i\ = 1): 

Re(y) > : Re [ Vl (y)] > , y e H : <p' t (y) > 

(A2) 

yeU+ : <pi(y)<£ln2 + 2y 

We assume that fi(y) behaves as a power-law for large y: 

y^oo : tpi(y) ~ a t y 1 ^ 1 with < W; < 1 (A3) 

The bounds on lui are due to (A2). In the following we 
only consider < u>i < 1 which is consistent with numer- 
ical simulations (Colombi et al.1997; Munshi et al.1999; 
Valageas et al. 1999a) which give 0.3 < u>i < 0.7 for 
— 2 < n < 1. Note that since we do not assume the coef- 
ficients S p (l) to be scale-invariant we write explicitely the 
I— dependence of the function ipi(y) and its parameters a/, 
loi . Then, we consider the case where the rapid growth of 
the coefficients S p with p implies a singular behaviour of 
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<pi(y) at small negative values of y, say y S: i = —l/x s j with 
x St i being large (x s j ~ 10): 



V^vt.i ■ Pi(y) = ~a s ,iT(uj S! i) (y - y S! iY 



(A4) 



where we neglected less singular terms. This leads to an 
exponential cutoff at large densities of Pi (p) which agrees 
with numerical simulations (Colombi et al.1997; Valageas 
et al. 1999a) and this singular behaviour of (fi(y) arises 
naturally from a tree-model for the correlation functions 
(Bernardeau & Schaeffer 1992). 
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Fig. 2. Probability Pi(p) for finding an overdensity p in a 
cell within which the average correlation function is £ = 200 
in the case uii = 1/2, u Si i = —3/2 and x St i — 10. The 
density distribution is a power-law with a cutoff at small 
(pv ~ £ "i/f 1 "i'-j an( j } ar g C [p c ^ densities. Both p v 
and p c are scale-dependent: for large £, p„ goes to while 
p c goes to infinity. The dashed line is the approximation (A6) : 
Pi(p)ocp^- 2 . 



This form of <pi(y) leads to a specific behaviour for 
Pi (p) in the non- linear regime £ 3> 1 (Balian & Schaeffer 
1989a). At very large overdensities one gets a pure expo- 
nential cutoff: 



p»l : P^^QTp" 'exp 



(A5) 



In the intermediate range Pi (p) is a power- law: 



£ ^ «p«C : P;(p) 



r(^)? 2 V? 



(A6) 



At very low densities ( "voids" ) one obtains an exponential 
cutoff: 



x cxp 



1 - Ull 



(A7) 



if 

where we defined z = p a t ^ ■ This be- 

haviour of the probability distribution Pi(p) is shown in 
Fig.2. 

Of course, since ipi (y) defines Pi (p) it also determines 
the moments p p (see Balian & Schaeffer 1989a, 1989b, for 
details). Using (21) one obtains: 



: {-y)-* My) e 



-vi(y)/Z 



(A8) 



with Re(y) < 0. For large p and £ one can drop the expo- 
nential: 



p> (1-wj) , 1 : 

/+ZOO 7 
loo 2^ ( ^ rP ^ (y) 

In particular, for p € IN* we recover (13): 

V e IN* , £ » 1 : Mp (Z) ~ S p (/) = ? P (0 

On the other hand, we can also write from (21): 



(A9) 



r- 1 



r(i-p) 7 

which leads to: 

p < 1 - , £ > 1 : 



(A10) 



(All) 



r l- 



(A12) 



'l J p/(l-u>i) J -pui/(l-ui) 



T(l-p) 



Note that (A12) gives the exact leading order in the limit 
p — > — oo at a fixed value of £, or in the limit £ — > oo at 
fixed p. 

References 

Balian R., Schaeffer R., 1989a, A&A 220, 1 
Balian R., Schaeffer R., 1989b, A&A 226, 373 
Bernardeau F., Schaeffer R., 1992, A&A 255, 1 
Bernardeau F., 1994a, A&A 291, 697 
Bernardeau F., 1994b, ApJ 427, 51 
Bertschinger E., 1985, ApJS 58, 1 



14 



P. Valageas: Non-linear gravitational clustering: smooth halos, substructures and scaling exponents. 



Bond J.R., Kofman L., Pogosyan D., 1996, Nature 380, 603 
Colombi S., Bouchet F.R., Schaeffer R., 1992, A&A 263, 1 
Colombi S., Bouchet F.R., Schaeffer R., 1994, A&A 281, 301 
Colombi S., Bouchet F.R., Schaeffer R., 1995, ApJS 96, 401 
Colombi S., Bouchet F.R., Hernquist L., 1996, ApJ 465, 14 
Colombi S., Bernardeau F., Bouchet F.R., Hernquist L., 1997, 

MNRAS 287, 241 
Davis M., Peebles P.J.E., 1977, ApJS 34, 425 
Frisch U., 1995, Turbulence, Cambridge University Press 
Ghigna S., Moore B., Governato F., Lake C, Quinn T., Stadel 

J., 1998, MNRAS 300, 146 
Hernquist L., 1990, ApJ 356, 359 

Jain B., Mo H.J., White S.D.M., 1995, MNRAS 276, L25 
Lacey C.G., Cole S.M., 1993, MNRAS 262, 627 
Lacey C.G., Cole S.M., 1994, MNRAS 271, 676 
McClelland J., Silk J., 1977, ApJ 217, 331 
Moore B., Governato F., Quinn T., Stadel J., Lake G., 1998, 
ApJL 499, L5 

Munshi D., Bernardeau F., Melott A.L., Schaeffer R., 1999, 

MNRAS 303, 433 
Navarro J.F., Frenk C.S., White S.D.M., 1996, ApJ 462, 563 
Navarro J.F., Frenk C.S., White S.D.M., 1997, ApJ 490, 493 
Peacock J.A., Dodds S.J., 1996, MNRAS 280, L19 
Peebles P.J.E., 1980, The large scale structure of the universe, 

Princeton University Press 
Press W., Schechter P., 1974, ApJ 187, 425 
Salvador-Sole E., Solanes J.M., Manrique A., 1998, ApJ 499, 

542 

Tormen G., Bouchet F.R., White S.D.M., 1997, MNRAS 286, 
865 

Valageas P., Schaeffer R., 1997, A&A 328, 435 
Valageas P., 1998, A&A 337, 655 

Valageas P., Schaeffer R., 1998, accepted by A&A, astro-ph 
9812213 

Valageas P., Lacey C, Schaeffer R., 1999a, submitted to MN- 
RAS, astro-ph 9902320 

Valageas P., Schaeffer R., Silk J., 1999b, accepted by A&A, 
astro-ph 9903388 

White S.D.M., 1979, MNRAS 186, 145 



This article was processed by the author using Springer- Verlag 
F/TeX A&A style file L-AA version 3. 



